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Abstract 

We study the phase transitions of the spin-± Ising film in a transverse field within 
the framework of the effective field theory. We evaluate the critical temperature of 
the film as a function of the exchange interactions, the transverse field and the film 
thickness. We find that, if the ratio of the surface exchange interactions to the bulk 
ones R = J s /J is smaller than a critical value R c , the critical temperature T c /J of 
the film is smaller than the bulk critical temperature T C B / J and approaches T® / J as 
R increases further. On the other hand, if R > R c , T c /J is larger than both the bulk 
I J and the surface T^/J critical temperatures of the corresponding semi-infinite 
system, and approaches T^/J as R increases further. 

P.A.C.S.: 75.70.Ak, 75.40.Cx, 05.50.+q 
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1 Introduction 



During the last years much effort has been directed towards the study of crit- 
ical phenomena in various magnetic layered structures, ultrathin films and 
superlattices [1-4]. The basic theoretical problem is the examination of the 
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magnetic excitation and the phase transitions in these systems. Of these, mag- 
netic films are very important both from the theoretical and the experimental 
standpoints [5,6], and can be studied as models of the magnetic size effect 
[7] and as quasi-two-dimensional systems. The magnetic and phase transition 
properties of semi-infinite Ising systems have been investigated for many years. 

The surface magnetism of these systems is very interesting [8-17]. It exhibits 
different types of phase transitions associated with the surface; if the ratio 
R, (R = J s /J), is greater than a critical value R c = (J s /J) C rit, the system 
may order on the surface before it orders in the bulk. As the temperature is 
lowered, the system undergoes two successive transitions, namely the surface 
and bulk phase transitions, whose critical temperatures are called the surface 
(T c s ) and the bulk (T C B ) critical temperatures respectively. On the other hand, 
if the ratio R is less than R c , the whole system becomes ordered at the bulk 
transition temperature T C B . 

Magnetic excitations in superlattices were considered in numerous papers (see, 
e.g., ref. [18] for a brief review). However, less attention has been paid to crit- 
ical behavior, and in particular to critical temperatures in superlattices. Ma 
and Tsai [19] have studied the variation with the modulation wavelength of 
the Curie temperature for a Heisenberg magnetic superlattice. Their results 
agree qualitatively with experiments on Cu/Ni films [20]. Superlattice struc- 
tures composed of alternating ferromagnetic and antiferromagnetic layers have 
been investigated by Hinchey and Mills [21,22], using a localized spin model. 
A sequence of spin-reorientation transitions are found to be different for su- 
perlattices with the antiferromagnetic component consisting of an even or odd 
number of spin layers. 

Fishman et al. [23] have discussed, within the framework of the Ginzburg- 
Landau formalism, the static and dynamical properties of a periodic multilayer 
system formed of two different ferromagnetic materials. They have computed 
the transition temperature and the spin-wave spectrum. On the other hand, 
the Landau formalism of Camley and Tilley [24] has been applied to calculate 
the critical temperature in the same system [25]. Compared to ref. [23], the 
formalism of ref. [24] appears to be more general because it allows for a wider 
range of boundary conditions and includes the sign of exchange coupling across 
the interface. 

For more complicated superlattices with arbitrary number of different layers in 
an elementary unit, Bernas [26] has derived some general dispersion equations 
for bulk and surface magnetic polaritons. These equations have been then 
applied to magnetostatic modes and to retarded wave propagation in the Voigt 
geometry [27]. 

On the other hand, with the developement of modern vacuum science and in 
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particular the epitaxial growth technique, it is possible to study experimen- 
tally the magnetic properties of low dimensional systems. For example, by 
depositing magnetic atoms on the top of non magnetic substrates, the thick- 
ness dependence of the critical temperature of ultrathin films of Gd on W(110) 
[28] and of Fe on Au(100) [29], has been measured. 

Our aim in this paper is to study the phase diagrams of an Ising (spin-±) film 
in a transverse field within the framework of the effective field theory [30]. 
This technique is believed to give more exact results than those of the standard 
mean-field approximation. In section 2 we outline the formalism and derive the 
equations that determine the layer magnetizations, the average magnetizations 
and the critical temperature of the film as functions of temperature, exchange 
interactions, transverse fields and film thickness. The phase diagrams of the 
film are discussed in section 3. The last section is devoted to a brief conclusion. 



2 Formalism 



We consider a spin-± Ising film of L layers on a simple cubic lattice with 
free surfaces parallel to the (001) plane, submitted to a transverse field. The 
Hamiltonian of the system is given by 

H = - E Ji^-^af, (1) 

<ij> i 



where of and of respectively denote the z and x components of a quantum spin 
<Ji of magnitude a — \ at site i, is the strength of the exchange interaction 
between the spins at nearest-neighbor sites % and j, and represents the 
transverse field acting on the spin at site i. We assume = J s if both spins 
belong to surface layers and Jij = J otherwise. 

The statistical properties of the system are studied using an effective field 
theory whose the starting point is the generalized, but approximate, Callen 
[32] relation derived by Sa Barreto et al. [33] for the transverse Ising model. 
The longitudinal and transverse magnetizations of the spin at any site % are 
approximately given by (for details see Sa Barreto and Fittipaldi [34]) 



1 



ml = (of) = - - 



tanh 




(2) 
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nf = {of) = -( r tanh f -(3[Q 2 + (V Jurf) 2 ]* 

= (fx (e Jitf, fii j ) = E ^) ) ; (3) 

where ml and are respectively the longitudinal and transverse magneti- 
zations at site i, (3 — l/k^T (we take k& = 1 for simplicity), (. . .) indicates 
the usual canonical ensemble thermal average for a given configuration, and 
the sum runs over all nearest neighbors of site i. We assume that the trans- 
verse field Q depends only on the layer index, which we shall denote by n. 
Because of the translational symmetry parallel to the (001) plane, also the 
magnetizations only depend on n. 

To perform thermal averaging on the right-hand side of eqs. (2,3), we follow the 
general approach described in ref . [30] . First of all, in the spirit of the effective 
field theory, multi-spin correlation functions are approximated by products of 
single spin averages. We then take advantage of the integral representation of 
the Dirac's delta distribution, in order to write eqs. (2,3) in the form 



% a n = J du f a (u, n n ) J dt exp(icut) JJJ (exp(-i£ J^*)) , (4) 
where a = z, x and 

W, An) = tanh + tfj) , (5) 

Uv, On) = ^-r tanh (i/% 2 + tfj) 

= f z (n n ,y). (6) 



We now introduce the probability distribution of the spin variables (for details 
see Saber [30] and Tucker et al. [31]): 



P{<) = \ 



(1 - 2m z n )5 (< + + (1 + 2m z n )5 (< - ±) 



(7) 



Using this expression and eq. (4), we obtain the following set of equations for 
the longitudinal layer magnetizations: 



N N 
fi=0 fj.i=0 
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x (1 -2m z 2 Y 1 (l + 2m z 2 ) N °->* 1 

x j z [i \R(N - 2fi) + (N - 2/iO] , fii) ; (8) 

N N N 

< = z- N - 2No E E E c»c»°c%(i - 2 ffl j)"(i + Kr" 

/i=0 /ji=0 /i2=0 

x (1 - 2r<_ i r (1 + 2r<_ 1 ) JV °-^(l - 2m* +i r (1 + 2< h i)"°-' la 
x / 2 [(TV + 2N ) - 2{fi + /i! + // 2 )] , , 

for n = 2,3,...,L-l; (9) 

mi = 2-^° E E t>C(! - + 2mi)^ 

/j.=0 £ti=0 

x (1 - 2m z L _ 1 Y 1 (l + 2m£_ 1 ) jV °-' il 

x / 2 [i [R(N - 2/i) + (N - 2^)] , ft L ) . (10) 

In these equations we have introduced the notation R = J s / J, N and iVo are 
the number of nearest neighbors in the plane and between adjacent planes 
respectively, and C l k are the binomial coefficients, C l k — l\/(k\(l — k)\). For the 
case of a simple cubic lattice, one has N = 4 and N — l. 

The equations for the transverse magnetization for each layer are obtained by 
substituting the function f x instead of f z in the expression of the longitudinal 
magnetization. Since however f x (y,Q) = f z (Q,y), this yields 

™>* = m z n [f g (y,n n ) — > f x (y, Sl n )] = m z n [f z (Q n , y)]. (11) 



We have thus obtained a set of self consistent equations (8-11) for the layer 
longitudinal and transverse magnetizations m z n1 m*, that can be directly solved 
by numerical iteration. However, since we are interested in the calculation of 
the longitudinal order near the critical temperature, the usual argument that 
the layer longitudinal magnetizations m z n should tend to zero as the temper- 
ature approaches its critical value, allows us to consider only terms linear in 
m z n) because higher order terms tend to zero faster than m z n . Consequently, all 
terms of order higher than linear in eqs. (8-11) can be neglected. This leads 
to the following system of equations: 

m z n = K,n-\m z n -i + A nn m z n + A njn+1 m z n+1 , (12) 



which can be written as 

K = <. (13) 
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where the matrix A is symmetric and tridiagonal with elements 

Aj = AAj + Aijtfij-x + 8 iJ+1 ). (14) 



For simplicity, we now assume that the transverse field acting on the system is 
uniform and equal to Q. In this case, the only nonzero elements of the matrix 
A are given by 



An = All = \ [f z (^(4R + 1), n) + /, (±(AR - 1), fl) 

+ 2f z [1{2R + 1), n) + 2f z (±{2R - 1), n) } ; (15) 

A 12 = A L)L _i = 1 {/, (i(AR + 1), n) - f z (±(AR - 1), n) 

+ Af z [i{2R + 1), n) - 4/, (|(2i2 - 1), n) + 6/, n) } ; (16) 

A nri = 4^ nin _! = 4A n ,„ +1 = 1 {/ 2 (3J, ft) + 4/ 2 (2 J, ft) + 5/ 2 ( J, ft)} ; 

for ra = 2,3,...,L-l. (17) 

The system of eqs. (12) is of the form 

M< = 0, (18) 



where the elements of the matrix M are given by 

Mij = (1 - Aii)5ij - AijiSij-! + 5 iJ+1 ). (19) 



All the information about the critical temperature of the system is contained in 
eq. (18). So far we have not assigned explicite values to the coupling constants 
and the transverse field: the terms in matrix (18) are general ones. 

In a general case, for arbitrary coupling constants, transverse field and film 
thickness, the evaluation of the critical temperature relies on numerical solu- 
tion of the system of linear equations (18). These equations can be satisfied 
by nonzero magnetization vectors m z n only if 

Det M = 0, (20) 
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where 



Det M = c 



a — 1 
-1 b -1 



-1 b -1 



-1 b -1 
-1 a 



(21) 



The parameters a, b and c that appear in equation (21) take into account the 
boundary conditions and represent the different propensities to order of the 
surfaces and of the bulk. They are given by 



I- A 



a = 



b = 



ii 



A 12 
1 — A 

A n ,n— 1 



l-A r 



L12 



l-Ar, 



A n ^n + 1 A nn /4: 

i \ 2(L-2) 



(22) 

for n = 2,3,...,L-l; (23) 



(24) 



In general, equation (20) can be satisfied for L different values of the critical 
temperature T c / J from which we choose the one corresponding to the highest 
possible transition temperature (cfr. the discussion in refs. [34,35]). This value 
of T c /J corresponds to a solution where ml, m^,. . . , m z L are all positive, which 
is compatible with a ferromagnetic longitudinal ordering. The other solutions 
correspond in principle to other types of ordering that usually do not occur 
here [34]. 

The reduction and rearrangement of the determinant of eq. (21) leads to the 
result [36-38] 



Det M = c[{ab - l) 2 D L _ 4 b - 2a{ab - l)D L _ 5 b + a 2 D L _ 6 b], 



(25) 
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where Dl(x) is the determinant 



D L (x) = 



x — 1 
-1 x -1 



-1 x -1 



-1 x -1 
-1 X 



(26) 



whose value is given by 



D L (x) = (x 2 -4)-5 2-( L+1 ) 



x 



X 



x- 



L+l 



X 



Vx^A 



L+l 



for x 2 > 4; (27) 



= sin[(L + l)k}/ sin k, with fc = cos _1 (a;/2), for x 2 < 4. (28) 

From now on, we take J as the unit of energy in our numerical calculations, 
and we measure length in units of the lattice constant. 



3 Phase diagrams 



From eqs. (20) and (25), we can obtain the phase diagrams of the film. The 
results show that there can be two phases, a film ferromagnetic phase (F), 
in which the longitudinal magnetization (m z = j Z)n=i m n) is different from 
zero, and a film paramagnetic phase (P), in which fn z = 0. In addition, if 
the number of layers in the film, L, is very large (L — > oo), the film should 
practically behave as a semi-infinite Ising system. It is well known that, if 
the ratio R = J s / J of surface to bulk exchange interactions in a semi-infinite 
Ising system, is greater than a critical value R c = (J s /J) C rit, there appear 
two different transitions: the surface transition and the bulk transition. The 
critical temperatures related to them are respectively called the surface critical 
temperature T c s and the bulk critical temperature T C B . To obtain the bulk and 
surface critical temperatures of the semi-infinite Ising system, we follow the 
approach due to Binder and Hohenberg [39]. Using eqs. (22-24), the system 
of linear equations (13) yields 



ami ~ m 2 = 0? 



(29) 
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—ml + bm z 2 - ml = 0; 
-m z n _ x + bm z n - m z n+l = 0, 



for n > 3. 



(30) 
(31) 



According to Binder and Hohenberg [39] , let us assume that m z +1 = 7 for 
n > 3, e.g., the layer magnetization m^ of each layer, with n larger than 2, 
decreases exponentially in the bulk. The equations (29) and (30) then yield 
the following secular equation: 



mi 



( 




0, 



(32) 



where, from eq. (31), the parameter 7 is given by 
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(33) 



Thus, the surface critical temperature T c 5 / J can be derived from the condition 
Det M s = 0, namely 



0(6-7) - 1 = 0. 



(34) 



We can now study numerically the physical properties of the surface and 
the bulk of the semi-infinite Ising system. Here it is worth noting that in 
our treatment the bulk transition temperature can be determined by letting 

m a = m n-l = m n+l = mZ * n e< i- (31), 1-e., 

6-2 = 0. (35) 



This yields 

f z (3J, n) + 4/2(2 j, n) + 5f z (j, n) = 8/3. (36) 

At I J = 0, eq. (36) yields the bulk critical transverse field value: f2f /J = 
2.3529. On the other hand, for the special case of the Ising model in the 
absence of the transverse field (ft — 0), eq. (36) reduces to 

tanh(3/9J) + 4 tanh(2/9J) + 5 tanh(/3J) = 16/3, (37) 

which is the Zernike [40] equation for the simple cubic lattice. The transition 
temperature is then determined as Tf/J = 1.2683. 

A useful expression for determining the critical value R c = (J s /J) cr it is there- 
fore given by the simultaneous solution of the equations (33) and (34). The 
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Fig. 1. The variation of R c = (J s /J) cr it as a function of the transverse field Q/J. 

variation of R c as a function of the transverse field is shown in Fig. 1. It shows 
that R c increases with the increase of the strength of the transverse field from 
its minimal value R™ m = 1.3069 for Q/J = and reaches its maximal value 
i^ nax = 1.3328 at the bulk critical transverse field Q^/J = 2.3529. 

We now calculate the (T c / J,R— J s / J) phase diagrams for different values of 
the transverse field and of the number of layers: typical results are shown in 
Fig. 2. Our phase diagrams are qualitatively different from the corresponding 
diagrams for the semi-infinite ferromagnet. The main difference is that we 
get in the film only one well defined critical temperature T c / J, instead of the 
two critical temperatures T C B / J and T c s /J. This temperature depends on the 
film thickness. According to these results we must give a new definition of 
R c . In semi-infinite systems R c was defined as the value of the parameter R 
above which the two critical temperatures T C B /J and T c s /J exist. However, 
according to Fig. 2 the parameter R c can now be defined as that particular 
value of R at which the critical temperature does not depend on film thickness 
(the cross-over point in Fig. 2). As it has been stated the numerical values of 
R c and related T c /J parameters are exactly the same as those found for the 
semi-infinite system. Furthermore, according to the definition of R c , it can 
be expected that the cross-over point in Fig. 2 should also define the critical 
temperature of the three-dimensional infinite bulk system, where the surfaces 
and the R parameter are of no importance. Fig. 2, where the bulk and the 
surface critical temperatures of the corresponding semi-infinite system are 
represented respectively by the dashed and dotted lines also shows that this 
is really the case. 
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Fig. 2. The phase diagram in the (T c /J, R = J s /J) plane, (a) VI/ J = 0; (b) VI/ J = 2. 
The dashed line is the bulk critical temperature T B /J and the dotted line is the 
surface critical temperature T c s /J of the corresponding semi- infinite Ising system. 

From Fig. 2a, which corresponds to the case when there is no transverse field 
acting on the system, we find that the value of R c corresponding to the cross- 
over point is equal to 1.3069, which is equal to the value reported by Wia- 
trowski et al. [42] and by Sarmento and Tucker [43]. For R < R c , the critical 
temperature T c / J of the film is smaller than the bulk critical temperature 
T B / J, and T c / J increases with the increase of L, approaching the bulk critical 
temperature T B / J = 1.2683 asymptotically as the number of layers becomes 
large. When R = R c , the critical temperature of the film T c /J is indepen- 
dent of L, and equal to T B /J. On the other hand, for R > R c the critical 
temperature of the film T c /J is larger than both the bulk T C B /J and the sur- 
face T c s / J critical temperatures of the corresponding semi-infinite system. The 
larger L, the lower T c / J, and, when the number of layers L becomes large, T B 
approaches asymptotically T c s / J. 

We now study the influence of the surface exchange interactions in the presence 
of a transverse field. Fig. 2b shows that a transverse field acting on the system 
increases the critical value of R c . For example, with a transverse field Q/J = 2, 
the critical value of R c is shifted from 1.3069 to 1.3199. At the same time of 
course for a given ratio of the exchange interactions R = J s /J, the critical 
temperatures of the film and of the corresponding semi-infinite system are 
reduced. 

In Figs. 3, we show the phase diagrams of the film in the (T c / J, VI/ J)-plane 
with different values of the thickness L, and the critical temperatures of the 
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Fig. 3. Phase diagram in the (T c /J, fi/J) plane: (a) R=l;(b)R= 1.5. The dashed 
line is the bulk critical temperature T B /J and the dotted line is the surface critical 
temperature T c s /J of the corresponding semi-infinite Ising system. 

semi-infinite Ising system for different values of the parameter R. The pres- 
ence of a transverse field, of course, reduces the critical temperatures of the 
film and of the semi-infinite system. We find that the (T c /J, Q/J) curve for 
a given value of R intercepts the (f2/J)-axis at a critical value Q c / J of the 
transverse field. When Q/J > Q c /J, there cannot be a ferromagnetic phase 
at any temperature. Fig. 3a, which corresponds to it! = 1 (smaller than R™ m ), 
shows that, for any finite value of the film thickness L, the film critical temper- 
ature T c / J is smaller than the bulk critical temperature T C B / J and increases 
as L increases, approaching T B / J for large values of L. Fig. 3b corresponds to 
R = 1.5 > R™ m , and shows that the film critical temperature T c /J is larger 
than both the bulk critical temperature T B / J and the surface critical tem- 
perature T c s / J of the corresponding semi-infinite system. T c j J decreases with 
the increase of L, approaching T c s / J for large values of L. 

We show in Fig. 4 the thickness dependence of the critical temperature of 
the film for different values of R and Q/J. Fig. 4a corresponds to the case 
when the transverse field vanishes. It shows that for any value of R below 
the critical value R c = i2™ m = 1.3069, the critical temperature of the film 
increases with L and approaches the bulk critical temperature T C B / J = 1.2683 
asymptotically as the number of layers becomes large. On the other hand, for 
R > R c , the critical temperature of the film becomes smaller as the number 
of layers increases, and approaches asymptotically, for large values of L, the 
surface critical temperature T c s /J which depends on R. In Fig. 4b, we show 
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Fig. 4. Thickness dependence of the critical temperature of the film for (a) VI/ J = 0, 
and (b) VI/ J = 1. The dashed line is the bulk critical temperature T C B / 'J and the dot- 
ted line is the surface critical temperature T c s /J of the corresponding semi- infinite 
Ising system. 

the thickness dependence of the critical temperature of the film when there is 
a transverse field acting on the system. We consider the case when VI/ J — 1 for 
several values of R. We see that Fig. 4b exhibits the same qualitative behavior 
as Fig. 4a, except that the transverse field reduces the critical temperatures. 

Finally Fig. 5 shows the variation of the critical transverse field Vt c / J as a 
function of the thickness of the film L for several values of R. The dashed and 
dotted lines correspond respectively to the bulk and surface critical transverse 
fields of the corresponding semi-infinite Ising system. For R < R c the critical 
transverse field of the film, Vt c / J, is smaller than the bulk critical transverse 
field Vl^/J, and increases with the increase of L, approaching Vlf/J for large 
values of L. For R > R c , Vl c /J is larger both than Vl^/J and Vl^/J, and 
decreases as L increases, approaching Vl^/ J for large values of L. 



4 Conclusion 



We have studied, within the effective field theory, the phase diagram of the 
transverse spin-± Ising film, where the exchange interactions between spins 
on the surfaces are different from those in the bulk. We have investigated 
the effects of the surface to bulk exchange interaction ratio, of the strength 
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Fig. 5. The variation of the critical transverse field of the film Q c /J at which the 
critical temperature of the film T c /J becomes zero as function of the number of 
layers L. The dashed line and the dotted line correspond respectively to the bulk 
and surface critical fields of the semi-infinite Ising system. 

of the transverse field and of the film thickness on the phase diagram. The 
results show that, in the film, there is only one critical temperature T c /J 
which depends on L, R and Q/J. We have identified a critical value R c of the 
parameter R, such that when R = R c , T c /J is independent of L; for R < R c , 
T c /J is smaller than T C B / J and for R > R c , T c /J is greater both than T C B /J 
and T c s /J. When L becomes very large, then for R < R c (R > R c ), T c /J 
approaches T C B / J (T c s / J) of the corresponding semi-infinite Ising system. 
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